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Radial propagation of geodesic acoustic modes 
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The GAM group velocity is estimated from the ratio of the radial free energy flux to the total free 
energy applying gyrokinetic and two-fluid theory. This method is much more robust than approaches 
that calculate the group velocity directly and can be generalized to include additional physics, e.g. 
magnetic geometry. The results are verified with the gyrokinetic code GYRO [J. Candy and R. E. 
Wahz, J. Comp. Phys. 186, pp. 545-581 (2003)], the two-fluid code NLET [K. Hallatschek and A. 
Zeiler, Physics of Plasmas 7, pp. 2554-2564 (2000)], and analytical calculations. GAM propagation 
must be kept in mind when discussing the windows of GAM activity observed experimentally and 
the match between linear theory and experimental GAM frequencies. 



I. INTRODUCTION 

Geodesic acoustic modes (GAMs) are axisymmetric 
poloidal ExB flows with finite frequencies coupled to up- 
down-antisymmctric pressure perturbations, which pro- 
vide the restoring force of the oscillation. Their frequency 
is proportional to ^/2cs/R with the sound speed defined 
as Cs = a/ (Tio -I- Teo)/mi. This can be shown analytically 
[l[ and has been confirmed in recent measurements in AS- 
DEX Upgrade 01 and Dlll-D [|. As GAMs are driven by 
the plasma turbulence, they are believed to play an im- 
portant role in limiting the turbulence strength and thus 
might give rise to stable equilibria with reduced energy 
transport [J. GAMs have attracted widespread atten- 
tion, coveriiig for example the generation and damping of 
zonal flows 0, @, 0] , the GAM spectrum (for circular flux 
surfaces) @ and GAM eigenmodes Q. In some works 
outward propagating GAMs are reported but no general 
rules determining the direction and speed of GAM propa- 
gation have been given 0, 0, [1] , a crucial point concerning 
the radial windows of GAM activity observed in ASDEX 
Upgrade @ or DIII-D Q and the match between linear 
theory and experimental GAM frequencies [Toj . 

The energy contained in a wave packet is transported 
with its group velocity as shown for example in Ref. pT| . 
Thus, one can calculate the group velocity of a GAM 
wave packet by compairing its total energy to its energy 
flux (Poynting flux). This provides a rather powerful 
tool for determining the direction and speed of GAM 
propagation, ft also allows a more general estimate of 
the maximal group velocity than a direct calculation of 
the dispersion relation, which requires much more effort 
and is restricted to a few simple cases. 

The basic concepts of the method are demonstrated for 
the two-fluid equations for cold ions and large safety fac- 
tor q. The generalization to warm ions, arbitrary safety 
factors and a gyrokinetic model is straightforward. The 
calculations are corroborated with exact analytical cal- 



culations for simple test cases and with numerical re- 
sults obtained by the gyrokinetic code GYRO and 
the two-fluid code NLET [l^] for various values of q and 
T = T,/T,. 

The structure of this article is as follows. In Sec. [Ill the 
equations for the free energy, its flux and the group veloc- 
ity are derived within the two-fluid framework for r = 0, 
q ^ oo, low j3, and high aspect ratio circular magnetic 
geometry. The calculation is generalized in Sees. Illll and 
IIVI to the warm ion, finite q case and to the gyrokinetic 
model. In Sec. |Vl the effects of plasma shaping on GAM 
propagation are discussed. Finally, the results are sum- 
marized in Sec. IVII 



II. 



FLUID MODEL FOR COLD IONS AND 
INFINITE SAFETY FACTOR 



The units are chosen such that the magnetic drift ve- 
locity is unity. Density, n, temperature, Ti and Tg, and 
electric potential perturbations (j) are normalized to 
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respectively, where the subscript indicates the corre- 
sponding background value and p* is given by Pse/R 
with the major torus radius R, Cse = (Tb.e/'Ti;)^/^, and 
Pse = (rriiCse) I (eS). The time scale is to = R/i'^Cse)- 

Beforehand, it is useful to recall the main character- 
istics of the GAM. An (m, n) = {Q,Q) E x B plasma 
flow, is generated by the flux-surface averaged electric 
field — V0O, with = {4>) and (...) indicating flux sur- 
face averaging. The divergence of the flow due to the 
magnetic inhomogeneities gives rise to a mainly up-down- 
antisymmetric TO = 1 pressure perturbation [Fig. [T](a)]. 
Since the compression of plasma requires work taken from 
the kinetic energy (V0o)^/2 of the flow, a restoring force 
is generated, the flow is slowed down, stopped, and even- 
tually reversed. An oscillation between pressure pertur- 
bations and flow results, in which the maximal energy 
stored in the pressure perturbations is comparable to the 
initial kinetic energy. 




(b) 
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FIG. 1: Sketch of a geodesic acoustic mode, (a) The poloidal E x B-flow, moving with the GAM phase velocity Vp, leads to 
compression or expansion of the plasma (indicated by the filled and striped areas, respectively) . Thus, an up-down antisymmet- 
ric m = 1 density perturbation arises, which is phase-delayed against the flow by n/2. (b) Density perturbations are associated 
with an E X B flow, which leads to compression or expansion of the plasma, causing a drift of the density perturbations drifting 
with the ion magnetic inhomogeneity drifts, (c) Parallel drift and phase velocity enhance the density perturbations caused by 
the GAM poloidal rotation, antiparallel Vd and Vp lower the perturbations. 



A. GAM Poynting flux and group velocity 

As a first approach, we calculate the Poynting flux (and 
group velocity) of the GAM for the cold ion two-fluid 
equations, neglecting sound waves (g oo). Since in the 
absence of perturbations the free energy is minimal, it 
is second order in the fluctuations. In the chosen frame- 
work, the free energy functional [l^ is given by 



{E) = {Ee + E,) = 




(V</)o)' 



(2) 



where E,, and Ei are the electron and the ion free energy 
density, respectively, v? /2 is the energy of the electron 
density perturbations, and (V0o)^/2 the ion kinetic en- 
ergy. 

The ion density fluctuations obey 



A(/) - C0 = 0, 



(3) 



where C = -Vd ■ V and Vd = -(1/2)(k + i?VlnS) x b 
is the sum of the curvature and Vi?-drifts of the elec- 
tron density fluctuations, and Ai/) is the divergence of 
the polarization current. The electrons are assumed to 
be adiabatic 



= ^ - 00, {n) = 0, 



(4) 



because the GAM frequency is much smaller than the 
electron bounce and transit frequencies. By combining 
([2]), ([3]), (HI), and representing the time derivative of {E) 
as the divergence of a radial Poynting flux one obtains 



dt (E)^~(W-S) = {-W- 



(5) 



The first term, Vdn,'^/2, represents the flow of the energy 
of electron pressure perturbations in ion magnetic drift 
direction. Since the radial component of Vd is up-down 
antisymmetric (for symmetric flux surfaces), this energy 
flux has a non-vanishing flux-surface average only if 
has an up-down asymmetry. That such an asymmetry 
exists and that the energy transport is in the ion drift 
direction can be shown as follows. 

Because due to adiabaticity ^ the pressure fluctua- 
tions shown in Fig. [T] (a) are connected to potential fluc- 
tuations, they are encircled by x B-flows as indicated 
in Fig. [1] (b). Similar to the poloidal flow, the vortices 
lead to compression or expansion of the plasma owing to 
the magnetic field variations. This effect is equivalent to 
the advection of the density perturbation by the ion cur- 
vature drift, computed with the electron temperature. 
Due to resonance between GAM phase and magnetic 
drift velocity, the pressure perturbations are enhanced at 
poloidal angles where drift and phase velocity are paral- 
lel, whereas they are weakened where those velocities are 
antiparallel. Therefore, an up-down asymmetry of the 
energy density arises, which leads to a net radial energy 
transport through the flux surface parallel to the phase 
velocity. Due to the asymmetry requirement, this flux 
somewhat resembles neoclassical density or temperature 
transport. 

Owing to adiabaticity pressure and potential per- 
turbations are equal, so that the gradient of the local 
density fluctuations causes an electric field, whose time 
dependence gives rise to a polarization current density 
— Vn. Since the temperature is normalized to Tg^ej the 
term — nVn can be interpreted as hydraulic energy flux 
Pejpoi consisting of the electron pressure Tq^^u and the 
polarization current density. We will refer to it in the fol- 
lowing as the polarization energy flux. As the E x _B-flow 



3 



associated with the density perturbations is proportional 
to the density gradients, the polarization energy flux can 
be regarded as the energy flux — ?;p(Vn)^, counterintu- 
itively with the reversed phase velocity Vp. 

The requirement of up-down asymmetry of the free 
energy density makes the curvature flux compara- 
ble in size to — nVri, a polarization effect. Whether the 
radial group and phase velocities eventually are parallel 
or antiparallel depends on the relative size of those two 
fluxes. 

Next, the free energy and the Poynting flux, Eqs. 
and (O, are evaluated in Fourier space for a circular high 
aspect ratio magnetic geometry. Recalling that the cur- 
vature operator C is up-down antisymmetric for circular 
flux surfaces, one obtains an estimate for the density per- 
turbations by splitting the ion density equation ([3]) into 
an up-down antisymmetric and a symmetric part with 
the corresponding densities Ua and 71^. Since the kinetic 
energy of the flow and the energy of the density fluctu- 
ations are of the same order, fc^^o ^ i only terms up 
first order in kr are kept in the antisymmetric equation. 
The symmetric density fluctuations are of second order 
in kr- Thus, with C = — sin(0)9,., the density becomes 



Vg,r 0.2 0.4 0.6 0.8 1.0 



ria w — sin (6*) (/)o. 



kr4>0- 



(6) 



Due to electron adiabaticity the GAM frequency u! deter- 
mined by (n) = (ris) = is given by uj = . Inserting 
© into ^ and ^ one obtains for the radial group ve- 
locity 



{E) 



(7) 



Since the ratio of curvature to polarization flux is —1/2, 
the total Poynting flux and the group velocity are an- 
tiparallel to the phase velocity for cold ions. 

The free energy approach only requires knowledge of 
the up-down antisymmetric density fluctuation Ua and 
its symmetric correction n^. The electron adiabaticity 
condition for n as given by Eq. ([6]) only yields the GAM 
frequency to 0th order in fc^. Thus, higher order cor- 
rections to the density have to be computed to calculate 
the group velocity directly from the dispersion relation. 
Hence, an advantage of the free energy approach is that 
less information is necessary compared to a direct calcu- 
lation of the GAM frequency. 

To verify the approximation ([7]), we give the exact so- 
lution of dl]), 



" [ sin (9) + n 



(8) 



with n = {Lj/kr){l + k^). The GAM frequency fol- 
lows from the condition {nsm(9)) = —Lijkr4>o, which 




FIG. 2: Exact GAM group velocity for r 
(solid) and approximation ((Tjl (dashed). 



m 



is obtained from Eq. ^ by using yielding oj = 
(2 -I- fc^)~^/^ and the corresponding radial group velocity 



(2- 



klf" 



2^2 



O [kl 



(9) 



which to lowest order in kr is identical to the approx- 
imation ([7]). Figure [2] shows the exact group velocity 
in comparison with the approximated result. For small 
wavenumbers the latter converges against the exact re- 
sult. Deviations for larger wavenumbers are due to drift 
velocity resonances. The resonance condition can be ob- 
tained from Eq. ([3]) giving 



UJ 



k ■ Vd 

1 + k^ 



sin (9) 
1 + fc2 



(10) 



for circular flux surfaces. A mode loses the character of 
a GAM, if, due to resonances, the energy of the den- 
sity perturbations becomes significantly larger than the 
kinetic energy. When the GAM frequency approaches 
a resonance, the density amplitude becomes very large 
compared to the poloidal rotation and dominates the 
mode. The resonant pressure perturbations propagate 
with the group velocity of the resonant drift mode. Ac- 
cordingly, ([3]) and ([6]) imply, that the modes discussed 
here are GAMs for small radial wavenumbers fc^ <C 1 
only. Therefore, the deviations of the approximate fre- 
quency from the exact one shown in Fig. [2] are due to 
the transition of the mode from a GAM to a magnetic 
drift mode. 



B. Numerical studies with NLET 

To corroborate the analytical insights, linear numerical 
studies were carried out with the two-fluid code NLET. 
The computations were performed on a grid of 1024 ra- 
dial and 32 parallel grid points with high aspect ratio 
circular geometry. The radial width of the computa- 
tional domain was 400 pse- The remaining parameters 
arc e = a/ R = 0, r = 0, and q ranging between 3 and 
30. Approximating e = a/R — 0, t = and q — > 00, 
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FIG. 3: NLET computed log-color-coded GAM spectrum 
(Fourier transform of 0) for t — and q = 30 with exact 
analytical frequency (solid), approximate frequency [dotted, 
obtained by integrating Eq. ((Tj] and magnetic drift resonance 
frequency kr/ (l -I- k^) (dashed). The exact frequency (bright 
line) interpolates between the approximate and the resonance 
frequency. Also visible are additional drift modes enclosed by 
the resonance line. 



the GAMs have been initiahzed at time t = with an 
(m, n) = (0, 0) electrostatic potential, which then evolves 
selfconsistently. 

The resulting spectral density of the radial E x B-How 
profile VEir,t) is shown in Fig. [3l The exact frequency, 
which also agrees with the numerical result, interpolates 
between the approximation obtained by integrating ([7]) 
and the resonance frequency. For kr < 0.5 the mode is 
obviously a GAM, whereas for larger kr it is gradually 
taken over by the magnetic drift resonance. For fc^ > 2 
the mode is completely dominated by resonant pressure 
perturbations as discussed in Sec. Ill Al and has lost the 
character of a GAM. 



III. GENERALIZED FLUID MODEL 

The generalization of the preceding calculations is 
straightforward. Following Ref. one can obtain the 
generalized ion free energy functional 



+ ^(VrT,f + ^[vl+r{Vv\\f), (11) 

which includes ion temperature and parallel velocity. For 
warm ions, the ion density and temperature fluctua- 
tions contribute to the total free energy with (l/2)Tn^ + 
(3/4)tJ;2. The FLR (finite Larmor radius) heat flux 
contributes the energy density (5/4) (VrT^)^. The en- 
ergy density of the diamagnetic drift velocity increases 
by rV(n -|- Ti). For finite safety factor, the parallel flow 
(l/2)w^ and the FLR correction (1/2)t(Vv||)^ also con- 



tribute to the total energy density. Fluid equations which 
exactly conserve the free energy functional are given 
by 



A 



1 - -tA dnv 



C (0 -f rn + rT, - rAa) = 0, 
2 
3 

C 



(12) 



■ Til + '^tTi 



1 - -tA 1 



• + m + 'i^TTi — rA5 



0, 



Aw'ii + - ^^a) (0 + r (n + T,)) 



0, 



(13) 



(14) 



with a = a{(j) + Tn) + PrTi and b = ^{ip + Tn) + jTTi. For 
collisionless plasma, the three coefficients a, /3 and 7 are 
given by (11/6,11/3,85/12). Inserting (fT2lfT4| into dtE 
and writing the result in terms of divergences, we obtain 
the Foynting flux 



dt{E) 



-V 



1 5 
Vd I -{n + Tp,f + -{TT,f + Tvl 



V 



{n + rpi) V (n + rpi) + TTi\74>Q 



3 



-(ac 



+ 2f3cd + jd^) -^[a (Vc)^ + 2/3VcVd 

+ 7 (Vd)') - (/)oAa + T (^Avf - 3 (V«||)')] }) , 

(15) 

in which c = (f) + rn and d = rTi. The first two di- 
vergences on the right hand side of Eq. (|15p represent 
the advection of the fluctuation energy by the magnetic 
drifts and the polarization drift in complete analogy to 
the first term on the right hand side of Eq. The last 
divergence in (fT51) is an FLR correction to the first one. 

We can approximate the two functionals pT|) . (|15p . and 
the group velocity by splitting the fluctuations in (|12II14[) 
according to their up-down symmetry and keeping only 
the lowest order terms as in Sec. Ill Al When the GAM 
frequency approaches the sound frequency, the ratio of 
the energy densities of parallel flow velocity and density 
perturbations to the ion kinetic energy densities increases 
and tends to inflnity close to resonance. Hence, the mode 
loses the character of a GAM. Sound wave resonance is 
negligible, if q ^ 1 (in practice q ^ 3 is sufficient). In- 
serting the approximate perturbations into pip and (jl5p . 
one obtains the group velocity 



2V6(3 + 5t)^/^ 

1 



[(-9 + 21t+189t2 + 265t3) 



4g2 (3 + 5r) 
-f5360r^ + 3675t'')] . 



135 + 1026r + 3324r2 



(16) 
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FIG. 4: Ion temperature dependence of the fluid [Eq. (US}] 
and the kinetic [Eq. (f23|) ] GAM group velocities for fcr = 0.1. 



(For details see appendix IB] ) All additional terms com- 
pared to Eq. ([7]) are positive, which causes the group 
velocity to change sign at r w 0.16 (Fig- S])- When 
calculating GAM eigenmodes as in Ref. [3], regions of 
evanescent and propagating GAM would be switched at 
this critical t because the group velocity is reversed. Ex- 
istence and properties of global GAM eigenmodes might 
be relevant for the efficiency of GAM excitation in the 
same way as they are for Alfven wave excitation [e.g., for 
toroidal Alfven eigenmodes (TAE)]. Since GAMs have 
recently been found to play an important role in nonlin- 
ear turbulence saturation [l5l |. their propagation might 
influence turbulent transport. 



IV. GYROKINETIC MODEL 

Generalizing the previous discussion to gyrokinetic 
theory we use the linear model [T^ 

a*/ + ^ • (r/ + Fo Jo0) + (r/ + FqJo^) = 



with the quasineutrality condition 



l-En 



- / Jo/d\- = 0. 



(17) 



(18) 



The velocity is the sum of the curvature and Vi? 
drift of the individual particles. Fq is the thermal back- 
ground distribution function, which is normalized such 
that /Fo d'^w = 1. Gyro-averaging is represented by the 
operator Jq and the thermal average of Jq is defined by 
To = / Fq Jq d'^u. The Fourier representations of Jq and 
Jq are Jo(t^/^wj^/s,.) and ro(fc,.) = exp{~Tkf.)Ia{Tkf) , re- 
spectively, with the Bessel function of the first kind Jq 
and the modified Bessel function of the first kind Iq . The 
ion free energy density is 



E, 



1 (rfY 



2Fn 



l-To 



(19) 



in which first term represents the energy of the fluctua- 
tions of the gyro-averaged distribution function, and the 



second one the energy of the gyrophase dependent fluc- 
tuations, i.e. the plasma polarization. Using (jl7p . dtE 
can be written as 



dt{E) 



T 2F, 



{^.VA-,n}^ 

{^•Vr/,n} +{r/,^.V0„} ^ 
It ) .j„ I t J Jo 

[nFo,^-V<^o}]d\,-l{<t>,^}^ 



(20) 



with K = rf + JquFo- The susceptibility operator 
X is defined by its Fourier transform x(fcr) = (1 ~ 
ro(A;r))/(rfc^). The brackets denote 



{a, b}j^ = a{K *b) — b {K * a) . 



(21) 



with * indicating convolutions. As shown in appendix 
1X1 {a, b}j^ can always be written as divergence, provided 
the kernel K is symmetric. 

The first term on the right hand side of Eq. (PP]) rep- 
resents the advection of the free energy of gyro-averaged 
fiuctuations by magnetic drifts. The remaining four 
terms in the integral are FLR corrections (e.g., gyrovis- 
cosity and the Bakshi-Linsker effect) equivalent to the 
FLR energy fluxes in Eq. The last two terms 

describe the polarization energy flux. They can be ex- 
pressed as 



Xn (uikr-n) — Vp 



X |fcr0| dh\x 
2 9 In /cr ' 



(22) 



where the first term is the gyrokinetic equivalent to the 
fiuid term — riVn, and the second one is an FLR correc- 
tion. Keeping only the lowest order terms, we compute 
an approximation of the group velocity as in Sees. Ill Al 
and mil by splitting / and n in and (j^O)) according 
to their up-down symmetry, and even and odd terms in 
. One obtains the radial group velocity (App. |B| 



[(-32 + 24t + 586t2 + mir^) 



8y2(4-h7r)^/^ 

^- ^ (640 -I- 8736r -f 55000t2 

4g2(4-H7r)^ ^ 

215268r3 + 560074t* + 526209t''^)] . 



(23) 



Far from drift and sound resonances, for fc^ ^ 1, g S> 
1 , the gyrokinetic equation can be solved alternatively by 
a power series expansion in terms of kj. and q yielding an 
identical result, also in agreement with [5]. The difference 
between the kinetic (j23|) and the fluid (|16p group velocity 
is negligible for g — ^ oo (Fig. 2]) and for t = 0. For g > 3 
and T > 0.2 the kinetic group velocity tends to be higher 
than the fluid one (75% at r = 1, increasing with r). This 
is caused by an earlier onset of the coupling to parallel 
modes due to hot particles. 
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FIG. 5: GYRO computed log-color-coded GAM spectrum 
(Fourier transform of 0) for t = 0.5 and large safety factor 
q = 30. The analytical kinetic frequency (|23|) is indicated by 
the dashed black line. Due to damping of the resonant modes, 
the simulated dispersion ends at kr ~ 0.4, as for warm ions 
the number of resonant particles increases with r, so that 
Landau damping has to be taken into account. 



Computational analyses with GYRO confirm the an- 
alytical results obtained from the kinetic and the fluid 
calculation. The simulations have been performed on a 
grid of 800 radial, 1 toroidal, 6 orbit and 8 energy grid- 
points with a radial box size of 400/9se, ^ t ^ 0.5, 
and 3 ^ (7 ^ 30. Beforehand, agreement with NLET for 
cold ions has been checked. An example for t = 0.5 and 
g = 30 is shown in Fig. O together with the frequency 
obtained by integrating Eq. 



V. MAGNETIC GEOMETRY EFFECTS 

Let us now turn to the effects of plasma shaping on the 
Poynting flux. Assuming, for simplicity, cold ions, infi- 
nite safety factor, and neglecting the polarization drift, 
one can approximate Eq. ^ a,t 9 = ±tt/2 (where kr, Vd, 
and Vp are parallel) 



2VE 



2VE 



{u ~ kr ■ Vd) R Rkr ■ {Vp - Vd) ' 



(24) 



in which R is the major radius, and ve = kr4>o is the 
E X B-drift velocity. For kr <^1, equivalent to Vp Vd, 
Eq. dill) implies 



4vi 



1 + 



2Vd.r 



(25) 



Vd,r being the radial component of Vd- Accordingly, the 
flux surface averages of the two terms on the right hand 
side of Eq. ^ (the Poynting flux) can be estimated by 



1 



Vp 



krV 



±{n'),kru{n^) 



kr 



(26) 



Since {E) w {v?) and I'd = 1 in the units defined in|TTl the 
group velocity is of order ©[/c^]- As the typical velocity 
scale of turbulent motion is Vdia, Vdia 3> Vd, and kr "C 1, 
GAMs generally propagate much slower than turbulence 
and the magnetic drifts. 

The specific magnetic geometry enters the calculation 
by means of the factors kr{9)/uj and kr{0)uj in the neo- 
classical and the polarization and FLR fluxes, respec- 
tively. Experimentally, only the radial wavenumber at 
the outboard midplane, kg, is known. However, the GAM 
pressure fluctuations, energy fluxes, and group velocity 
have to be estimated at 6 = ±7r/2, where kr is smaller 
due to, e.g., ellipticity or Shafranov shift. For an el- 
liptic Miller equilibrium [l3|, the radial wavenumber at 
6* = 7r/2 is given by 



1 + drR 
K + rdrK 



kn 



(27) 



where k is the elongation and r the minor radius at the 
outboard midplane (r and R refer to the flux surface 
center). 

Typical values of the geometry parameters are [l7j 
drK = (k — l)/r, drR = —1/3 and aspect ratio A = 3.5. 
The K, dependence of lo for fc^ = can be obtained numer- 
ically [l^l and is rather accurately described for k = 1..2 
and g > 1 by 



LO (k) 



1 

V 



5 2.97 
3^1 + 3.73k' 



(28) 



Other possible parametrizations are discussed in Ref. 0] ■ 
Substitution of these parameters into Eq. ([?7|) yields 



3 (2k- 1) 



ko- 



(29) 



Multiplying the magnetic drift energy fluxes (|BmB3p 
with krUj{l) / {kQUj{K)) and the polarization energy fluxes 
(|B4IIB8|) with kfuj{K)/{k'^uj{l)), and defining k=1 + 6k, 
we obtain an approximation of the group velocities at the 
outboard midplane for t = 



{5k) 



ko 



0.32 + 2.155k + 5.31(5k2 + 5.5A5k^ - 
[-0.050 - 0.0215k + 0.195k^ + Q.0626k^ 
0.062 + 0.265k + 0.335k2 + O.II^k^ 



25k4 



(30) 



and T = 1 



Vg^r (<5k) 



ko 



0.32 + 2.155k + 5.315k2 + 5.545k3 + 25k'^ 
[0.59 + 1.845k + 1.575k2 + 0.525k3 



0.54 + 2.515k + 3.455k2 + 1.145/. 



(31) 



Figure [S] shows the k dependence of the group velocity. In 
the cold ion case, increasing plasma elongation leads to 
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6: Estimates ol the radial group velocities for kr — 
r = (solid) and r = 1 (dashed) plotted against the 



a change of sign of the group velocity when the "neoclas- 
sical" fluxes become larger than the polarization terms. 
Compared to circular flux surfaces, the group velocity 
for warm ions is reduced in elliptic geometry. There- 
fore, far from resonances, Vg ~ krVd remains small com- 
pared to the diamagnetic velocity and the magnetic drift. 
Our cold ion approximation agrees with numerical stud- 
ies performed with NLET. 

However, in single-null configuration near the separa- 
trix, the perturbations vanish at the X-point due to the 
magnetic null and GAMs are located opposite towards 
the X-point. Consequently, the neoclassical energy flux 
is (wd,r"-^)/2 ~ Vd(n^) and Vg ~ Vd, because the polariza- 
tion energy flux is one order smaller and can be neglected. 
Hence, independent of the ion temperature, GAMs prop- 
agate in the ion magnetic drift direction, which is usually 
directed towards the X-point, i.e. radially inward. The 
GAM dispersion must be linear in this case. In the up- 
down symmetric case, the radial mode structure can show 
standing waves while in up-down antisymmetric geome- 
try propagating waves are expected as the wavenumbers 
of incoming and reflected wave at a cutoff layer are dif- 
ferent. 



VI. SUMMARY AND DISCUSSION 

We have given estimates for the radial group veloc- 
ity of GAMs using an energy approach, which exploits 



Poynting's theorem. 

In the first step, the GAM energy flux ([5]) has been cal- 
culated in a two-fluid framework for cold ions and inflnite 
safety factor, and two kinds of transports have been iden- 
tifled. The "neoclassical" energy flux, which represents 
the advection of the free energy of the fluctuations with 
the magnetic inhomogeneity drift, {vd,rEfiuc) , is always 
parallel to the phase velocity. However, the polarization 
energy flux, {—k^pl^VpEjiuc), is always antiparallel to 
the phase velocity. Therefore, the one of the two fluxes 
which is dominant controls the propagation direction of 
the GAMs. This result has been generalized to include 
ion temperature and parallel flows (|15|) and to a gyroki- 
netic model (j20[) . in which equivalent energy transports 
have been found. 

In the second step, the group velocity, which is given 
by the ratio of Poynting flux to total free energy, has 
been evaluated for circular high aspect ratio flux sur- 
faces [Eqs. ([T5)) . and estimated for elliptic plasmas 
[Eqs. dSOl), JSH)]. The results agree with NLET and 
GYRO computations, and analytical calculations. For 
cold ions, group and phase velocity have been found to 
be opposite. For t > 0.2 or k > 1.5, they are paral- 
lel. The neoclassical energy flux requires an up-down 
asymmetry of Efi^c to give non-zero values. In case of 
up-down symmetric flux surfaces, the asymmetry is of 
order krPseEfiuc, which makes the neoclassical compara- 
ble to the polarization energy flux. Therefore, the group 
velocity Vg^r = {Sr)/{Efiuc) is of order krPseVd- 

Due to resonances with magnetic drift modes (not drift 
waves) and sound waves, the propagation speed of GAMs 
is limited. Close to such a resonance, the poloidal rota- 
tion becomes negligible compared to the remaining de- 
grees of freedom, as the amplitudes of density and tem- 
perature, or parallel flow velocity diverge. The charac- 
teristics of the resulting mode are completely determined 
by the resonant mode, and it should not be called GAM 
any longer. The resonances restrict GAMs to fc^ ^ 1 and 
1 and limit the group velocity. 

Since for single-null configuration, which is quite com- 
mon in today's experiments, the fluctuations vanish at 
the X-point, the neoclassical energy flux is of order 
VdEfiuc, and GAMs propagate at the magnetic drift ve- 
locity. More precisely, keeping only the first term on 
the right hand side of neglecting v\\, and using 

T ~ 2n/3, the ratio of Poynting flux to free energy yields 
the group velocities Vg^r ~ Vd for r = and Vg_r ~ Sw^ for 
T = 1. Since usually Vd is directed towards the X-point, 
GAMs propagate radially inward in this case. Overall, 
the GAM group velocity has been shown to be much 
smaller than the diamagnetic velocity, which is the typ- 
ical scale of turbulent motion. In up-down symmetric 
magnetic geometries, Vg^r is even much smaller than the 
magnetic inhomogeneity drift. 

Calculating the group velocity by means of the Poynt- 
ing flux is advantageous compared to a direct calculation 
of the GAM dispersion relation because only lowest or- 
der approximations of the fluctuations are required (Sec. 
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III A[) to gain insights into effects induced by ion tem- 
perature, sound waves, or magnetic geometry. Requiring 
relatively small efforts, the energy approach enabled us to 
calculate the group velocity rather accurately for circu- 
lar high aspect ratio geometry, to estimate it for elliptic 
Miller equilibria and to predict the effect of the X-point in 
single-null configuration. Since the formation and qual- 
ity of the H-modc is dependent on whether the magnetic 
drift is directed towards or away from the X-point, one 
may speculate if there is a relation to GAM propagation. 



iVQii + brf'^ [3 + 20r 
(81 + 387t + 500t2) 
^ V (3 + 5r) 

_ rfc302 (9 + + 50r2) 
(72^6 (3 + 5t)^/' 



, (B2) 



(B3) 



APPENDIX A: GYROKINETIC COMMUTATORS 
AS DIVERGENCES 

The brackets defined in Eq. (|2ip in Sec. llVl can always 
be written as a divergence. Consider the expression 



C ^a{K*h)-b{K*a) 



(Al) 



{ac^ + 2(3cd + ^(f) + 00 Aa 

+ (Vc)^ + 2/3VcVd + 7 (Vd)^) 

- T (^Ai;| - 3 (Vi^ii' 

- 1 \ llrfc302(3 + 7^) 

V; 3V6(3 + 5r)^/2 ' ^ ' 



where a{x) and b{x), K [x) ~ K {—x) is a symmetric con- 
volution kernel and * indicates the convolution operation. 
Because of this symmetry one may express equation (|Aip 



(- {n + Tpi) V (n + Tpi)) 



= - I 



1 ^^ fc302 (3 + 5^)3/2 

3V6 



, (B5) 



d,J = C = d, 



L {z — X, z — y) a{x)b (y) dx dy 



(A2) 

with L {x, y) = [9 (x) — 9 [y)] K (x — y). For purely har- 
monic waves a{x),b{x) = {(A, S) cxp (i/cx — ia;^)}, 
the resulting flux J is given by 



J = [d,)-^ {a, 6}^ = 3 {A*B} dkk (k) 



(A3) 



where K is the Fourier transform of K. The application 
of this theorem to, for instance, the first bracket in Eq. 
(I20[) yields the corresponding free energy flux 



Jflr 



Ji (fc±) ^^5R {n* {tJ + J^o Jon)} d\. 



(A4) 



APPENDIX B: INDIVIDUAL POYNTING 
FLUXES 

Evaluating Eq. (|15p for circular flux surfaces by using 
estimates of the up-down symmetric and antisymmetric 
fluctuation amplitudes (as described in Sec. Ill A[) the 
individual Poynting fluxes result up to order 0[fc3], 



1 + 1°-^^;^%,, , (B6) 



-tT,V0o = - 



3^/6(3 + 5T)^ 



V/ 3V6(3 + 5t)'/' 



(B7) 



For the gyrokinctic framework, Eq. ((20)) . the individ- 
ual fluxes arc 



T 2F, 



d^w ) = 







481r2 + 747r^ 



4\/2(4 + 7t)^/' 
1 



g2 (4 + 7r)- 



[16-f 140t 



(576 



6576r + 36068t2 -f 130144r3 

+ 317687r'* + 293067r^)] , (B9) 



(-'(-/{v--."},/S) 



0, (BIO) 



— [n + Tp,) 



Z-^" [9 + 30t + 55t^ 
6^/6 (3 + 5t) ^ 

^ (81 + 342t + 455t2) , (Bl) 



^/2Tfc302 (2 + 5r) 



g2 (4 7r) 



3/2 



(Bll) 
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2V2(4 + 7r) 



1/2 



13 



28 + 196r + 163r^ 



g2 (4 + TtY 



(B12) 



2V2 (4 + 7r) 



1/2 



12 



48 + 168t + 60t^ 
92 (4 + Irf 



(B13) 



-16 + 12r + 21r^ 



8^/2(4 + 7r)'/' [ 92(4 + 7^)^ 

(-192 + 656r - 252t2 + 612r3 + 483t'^)] . (B14) 
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